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Abstract 

Assuming that the matter filling the background geometry in the 
Early Universe was a free gas and no phase transitions took place, we 
discuss the thermodynamics of this closed system using classical ap- 
proaches. We found that essential cosmological quantities, such as the 
Hubble parameter H, the scaling factor a and the curvature parameter 
k, can be derived from this simple model. The results are compatible 
with the Friedmann-Robertson- Walker model and Einstein field equa- 
tions. Including finite bulk viscosity coefficient leads to important 
changes in the cosmological quantities. Accordingly, our picture about 
evolution of the Universe and its astrophysical consequences seems to 
be a subject of radical revision. We found that k strongly depends 
on thermodynamics of the cosmic background matter. The time scale, 
at which negative curvature might take place, depends on the relation 
between the matter content and the total energy. Using quantum and 
statistical approaches, we introduced expressions for H and the bulk 
viscosity coefficient. 

1 Introduction 



We assume that the cosmological background geometry filled with matter 
or radiation is characterized by well known physical laws, like equation of 
states and thermodynamics. It leads to a solid physical description for the 
Universe in its very early stages. Such a description, especially in the Early 
Universe, is favored, because we so far have no observational evidence against 
it. Other components of the cosmological geometry, like dark matter and 
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dark energy wouldn't matter much during these early stages. Therefore, we 
can disregard them. 

In this work, we introduce a toy model based on thermodynamical ap- 
proaches to describe the Early Universe. We disregard all phase transi- 
tions and assume that the matter filling the background geometry was likely 
formed as free gas. We apply the laws of thermodynamics and fundamentals 
of classical physics to derive expressions for the basic cosmological quanti- 
ties, such as the Hubble parameter H(t), the scaling factor a(t) and the 
curvature parameter k. We compare them with the Friedmann-Robertson- 
Walker (FRW) model and Einstein field equations. 

In this treatment, we apply the standard cosmological model and use 
natural units in order to gain global evidences supporting the FRW model, 
although we disregard the relativistic and microscopic effects. The var- 
ious forms of matter and radiation are homogeneously and isotropically 
distributed. We use non-relativistic arguments to give expressions for the 
thermodynamic quantities in the Early Universe, which obviously reproduce 
essential parts of FRW model. We assume that the Universe was thermal 
equilibrium and therefore the interaction rates exceed the Universe expan- 
sion rate, which was slowing down with the time t. Also, we assume that 
the expansion was adiabatic, i.e., no entropy and heat change took place. 

Finally, we take into consideration two forms of the cosmic background 
matter. The first one is ideal gaseous fluid, which is characterized by lack of 
interactions and constant internal energy. The second one is viscous fluid, 
which is characterized by long range correlations and velocity gradient along 
the scaling factor a(t). 

2 Expansion Rate in Non- Viscous Cosmology 

We assume that all types of energies in the Early Universe are heat, Q. In 
such a closed system, the total energy is conserved, i.e., 

dQ = = dU+pdV, (1) 

where U is internal energy, p is pressure and V stands for the volume. V 
can be approximated as a cube with sides equal to the scaling factor a, i.e., 
V = a 3 or as a sphere with radius equal to a, i.e., V = (47r/3)a 3 . In both 
cases, V oc a 3 . Apparently, Eq. ((TJ is the first law of thermodynamics. The 
expansion of the Universe causes a change in energy density p = U/V, i.e., 
surely decreasing, which can be given as dp = dU /V — UdV/V 2 . In comoving 
coordinates, U is equivalent to the mass m and consequently to the energy. 
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From Eq. ([I]), we get 



dp = -3(p + p) — . (2) 
a 

Dividing both sides by an infinitesimal time element dt results in 

p = -3(p + p)H, (3) 

which is the equation of motion from FRW model, which strongly depends 
on the thermodynamic quantities, p and p, i.e., the equation of state (EoS). 
One dot means first derivative with respect to the time t. H is the Hubble 
parameter which relates velocities with distances; H = a/ a. 

The radiation-dominated phase is characterized by p = p/3 and therefore 
Eq. ([3]) leads to p oc a~ 4 . In the matter-dominated phase, p « p and 
therefore p oc a~ 3 , i.e., p oc V^ 1 . The energy density p is a function of 
temperature T. Then, we can re-phrase the proportionality in the radiation- 
dominated phase as p oc TV . 

Neglecting both cosmological constant A and curvature parameter k, 
we simply get that H 2 oc p. Then, the scaling factor in the radiation- 
dominanted phase a oc t 1 ^ 2 and in matter-dominanted phase a oc i 2 / 3 . The 
results are depicted in Fig. [TJ 



3 Expansion Rate in Bulk Viscous Cosmology 

Let us assume that one particle of mass m is located at a distance a from 
some point in the Universe. Such a particle will have, in the radial direction, 
kinetic energy ma 2 /2. In the opposite direction, it is affected by a gravita- 
tional force due to its mass m and the mass inside the sphere M = (47r/3)a 3 p. 
Then the particles's gravitational potential energy is —GMm/a, where G is 
the Newtonian gravitational constant. The total energy is 

E = l -rad 2 -G^, (4) 
2 a 

which can be re-written as 

a z + k = —Gpa 2 , (5) 
3 

Last equation is, exactly, the Friedman's first equation with curvature pa- 
rameter k = —2E/m, which apparently refers to negative curvatures as- 
sociated with various geometrical forms depending on both total energy E 
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and mass m. In the Friedman's solution, k can be vanishing or +1 or — 1, 
referring to flat or positively or negatively curved Universe, respectively pQ. 
Our toy model agrees well with the negatively curved Friedman's solution, 
especially when the particle mass m equals two times the total energy, i.e., 
m = 2E. 

According to recent heavy-ion collision experiments [2] and lattice QCD 
simulations [3] , the matter under extreme conditions (very high temperature 
and/or pressure) seems not to be, as we used to assume, an ideal free gas. It 
is likely fluid, i.e., strongly correlated matter with finite heat conductivity 
and viscosity coefficients (bulk and shear) [4j. Therefore, it is in demand to 
apply this assumption on the background geometry in Early Universe. The 
cosmic background should not necessarily be filled with an ideal free gas. In 
previous works [UE117], we introduced models, in which we included finite 
viscosity coefficient. The analytical solution of such models is a non-trivial 
one O El [7] . In the present work, we try to approach the viscous cosmology 
using simple models, in which we just use classical approaches. As we have 
seen, the classical approaches work perfectly in the non-viscous fluid. It is 
in order now to check the influence of viscous fluid on the cosmological evo- 
lution. The simplicity of these approaches doesn't sharpen the validity of 
their results. Surely, it helps to come up with ideas on reality of the viscous 
cosmology. 

We now assume that the test particle is positioned in a viscous surround- 
ing. Then the total energy, Eq. (|4|), gets an additional contribution from 
the viscosity work, which apparently resists the Universe expansion, 

E = -ma 2 -G rja 3 -, 6 

2 a a 

where r] is the bulk viscosity coefficient. We assume that the expansion of 
the Universe is isotropic, i.e., symmetric in all directions. Consequently, the 
shear viscosity coefficient likely vanishes. Comparing Eq. ([6]) with the Fried- 
mann's solution leads to another expression for the curvature parameter, 

* = _*>*«.. (7) 

m ma 

Comparing Eq. (JT]) with the three values of k given in the FRW model 
{k = +1,0, —1) [1], results in three expressions for the expansion rate a in 
the bulk viscous cosmology. 
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When k = +1, then the expansion rate or velocity reads 

{2E + m\ 1 . x 

a = (8) 

\ rj J a z 

It is positive everywhere and inversely proportional to a 2 . It doesn't depend 
on the comoving time t, directly. Apparently, its t-dependence is embedded 
in the f-dependency of E, m and rj. The scaling factor itself is 

2E + m\^ 3 A/3 



a(t) = (3— (9) 

In Fig. [TJ we compare this result with the non- viscous fluid as given in sec- 
tion [21 i.e., a(t) oc i 1 / 2 for radiation-dominated phase and a(t) oc i 2 / 3 for 
matter-dominated phase. For simplicity, we assume that all proportionality 
coefficients are equal. It is clear that the scaling factor in the viscous cos- 
mology is the slowest one. This would refer to the fact that the viscosity 
likely resists the Universe expansion. Increasing rj shrinks or shortens a(t), 
Eq. ([9]). At very small t, the expansion of the bulk viscous Universe is much 
rapid than the other two cases (non- viscous) . 

From Eqs. © and Q, the Hubble parameter reads 

Hit) = i. (10) 

Apparently, H doesn't depend on any of the thermodynamic quantities. It 
is always positive and decays with increasing t. 

When k = — 1, the expansion rate or velocity takes the form 

IE -ml 
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rj a 2 

It is only positive, i.e., the Universe is only expanding as long as m < 2E. 
Otherwise, the expansion rate or velocity decreases. Eq. (IllD sets the limit 
of the Universe contraction. This limit is reached, when the mass m exceeds 
twice the total energy E. In Early Universe, E » m and, consequently, the 
Universe started explosively, although k could have a negative value. Much 
later, E decreases according to this expansion and the matter production, 
meanwhile the mass m gains more and more contributions. At a certain 
point, the expansion rate turns to the backward direction. It is necessarily 
to mention here that, this toy model takes into account the visible energy 
and matter components only. The invisible components are not included in 
it. 



5 



3 



matter 

radiation - - - 
viscous 




Fig. 1: Scaling factor a as a function of comoving time t is depicted for vis- 
cous fluid (solid), radiation-dominanted (long dashed) and matter-dominanted 
(dotted) phases. 



The scaling factor also depends on the total energy E and mass m, 

a(t) = ^3 J t 1/3 . (12) 

a is positive as long as m < 2E. Otherwise it switches to negative values. Its 
time dependence, a oc i 1 / 3 , looks like the previous case at k = +1, Eq. ([9|). 
From Eqs. (llip and (112|) . the Hubble parameter reads 

H») = i. (13) 

As in the previous case, k = +1, H is always positive and doesn't depend 
on any of the thermodynamic quantities. 

When k = 0, the expansion rate or velocity will be 

IE 1 

a = ; 

2 r\ a 



(14) 



In flat Universe, a does not depend on m. It increases with increasing the 
total energy E and decreasing the viscosity coefficient rj. Also, the scaling 
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Fig. 2: The Hubble parameter H as a function of comoving time t is depicted 
for viscous fluid (solid) and ideal gas. The radiation- and matter-dominanted 
ideal gasses are drawn as long-dashed and dotted lines, respectively. 



3£V /3 + i/3 



rate, 



5-j f", d5) 

depends on E and rj, only. It doesn't depend on the mass m, i.e., the mass 
production doesn't affect the scale factor or the expansion. From Eqs. (TH1) 
and (USD 

*w = j- t - do) 

In Fig. [2j we depict from this model and compare it with the two 
cases when the background matter is a non-viscous gaseous fluid. The latter 
is likely dominated by radiation or matter, where H = l/2t and H = 2/3t, 
respectively. We notice that H in the viscous cosmology is faster than H in 
the non-viscous cosmology. Therefore, we conclude that the bulk viscosity 
causes slowing down the Universe expansion. 
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At k = 0, the total energy is likely dominated by the viscosity work 



E = -rja 3 -. (17) 
a 



Plugging this equation into the scaling factor, Eq. (|T5|) . leads to deceleration 

2d 

4 = ST < 18 > 

Two dots refer to second derivative with respect to time t. 



So far, we conclude that filling the cosmic background geometry with 
bulk viscous fluid strongly affects (moderates) the expansion rate of the 
Universe. The evolution of the scaling factor a is damped, when rj increases. 
The curvature parameter k, which appears as a constant in FRW model and 
Einstein equations p], depends on the total energy E, the particle mass m 
and the viscosity coefficient rj. In the positively curved Universe, a increases 
with increasing E and m. In the negatively curved Universe, a increases 
only as long as 2E > m. Otherwise, it decreases causing Universe contrac- 
tion. We conclude that the scaling factor in the flat Universe depends on 
E/rj, but not on the mass constent, m. 



4 Energy Density in Bulk Viscous Cosmology 

Including the work of bulk viscosity into the first law of thermodynamics, 
Eq. ([1]), results in 



dU = - \pdV + r)a a -rdV) . (19) 

Apparently, the evolution of energy density depends on the Hubble param- 
eter H and the viscosity coefficient r), 

3(p + p) + 3r]a 3 ^j H. (20) 

Comparing this evolution equation with the one in the Eckart [8] relativistic 
cosmic fluid leads to a direct estimation for the bulk viscous stress IT. The 
conservation of total energy density requires that the bulk viscous stress 
equals to the work of bulk viscosity, i.e., II = rja 3 dd/da. 
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In the radiation-dominanted phase, EoS reads p = p/3 and Eq. (f20|) 
can be solved in comoving time t by utilizing our previous result on H(t), 
Eq. (fT6|) for instance, 

p = -lnt^ 12p ^ 3 l (21) 

On one hand, it implies that p diverges at t = 0. On the other hand, the 
viscosity coefficient r\ seems to moderate the evolution of the total energy 
density. Obviously, this result strongly depends on EoS [5j [7] , which is 
different in the different phases of Early Universe, i.e. differs with t. 

Eq. (|20p is consistent with the second law of thermodynamics at non- 
negative entropy production, S\ = II 2 jrfT > 0. In this model, S l i oc VH/T. 
The Friedmann's second solution in flat Universe, whose background geom- 
etry is filled with a non- viscous fluid, 

n 4-7T 

- = -^-G{p + p). (22) 
a 3 

seems to follow the second law of thermodynamics. To show this, let us 
take the time derivative of last expression. Then d(pa s ) = —a 2 (p + p)da. 
Depending on EoS, for instance in de Sitter Universe, last equivalence can 
be re-written as 

d(pa 3 ) = -3a 2 pda. (23) 

It is nothing but the second law of thermodynamics {dU = —pdV + TdS) of 
an adiabatic system, i.e., the expansion is thermally reversible and obviously 
doesn't affect the entropy content, dS = 0. 



5 Hubble Parameter in Quantum Cosmology 

Let us suppose that N particles are adhered to a cubic or spherical volume, 
i.e., V oc a 3 . The particles are distributed according an occupation function, 
which depends on their quantum numbers and correlations. According to the 
standard cosmological model, they are allowed to expand in a homogeneous 
and isotropic way. We suppose that particles have no interactions. Then, the 
energy of a single particle E = (p 2 -\-m 2 ) 1 / 2 , where momentum p = -(nix + 
n2y + ri3z). In natural units, h = 2ir. The state density in momentum space 
a 3 //i 3 = V/(2tt) 3 . From the integral of particle density in phase space, we 
get the particle density in ordinary space and, therefore, 

o 2vr 2 / f°° p 2 dp . 

« 3 = N — / / /_/ , (24) 
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where [i is the chemical potential and g is the degeneracy factor. Taking the 
time derivative (equivalent to 1/T) and dividing both sides by the scaling 
factor a results in 

1 fl g f°° 11 -E 2 , \ dT 



H ™ = ^{iwj. i^v/ */*- <25) 

where n is the particle number density, n depends on the intensive state 
variables, T, p and \i. It implies that H depends on 1/nT 2 and the time 
derivative, dT/dt, besides the integral, which can be calculated, numerically, 
in dependence on T and [i. It is obvious that the expansion of the Universe 
is driven by generating new states. 

When assuming that the background geometry is filled with a relativis- 
tic Boltzmann's gas, then the equilibrium pressure and energy density at 
vanishing viscosity are given as 

p(m, T) nT, 

/3T K 1 (m/T)\ . 
em,T « n — + — — =r m, 26 
V m K2{m/1 ) J 

where is the i-th order modified Bessel function. At equilibrium, the en- 
tropy is maximum. At vanishing chemical potential, the Hubble parameter 
in Eq. (j25l) reads 



H = 9 L 

2tt 2 6n 



-1 ( tyT \ (2 . orr 2\ 



pE + 4T.Mtanli -f— - (m z + 8T Z ) ln(2p + 2E) 
V EA4 J 

where A4 = (m 2 + 8T 2 ) 1//2 . In the relativistic limit, i.e., m — > 0, 
9 1 



dT 
dt ' 



H 



2vr 2 6n 



2 + 8T 2 ( v^tanh- 1 ^ - ln(4p))l ^. (27) 



We conclude that H in non-viscous quantum cosmology H depends on the 
intensive state quantity T, its decay with the time t, state density in mo- 
mentum space and both of momentum and number of occupied states [S]. 



6 Bulk Viscosity in Quantum System 

In this section, we give estimates for the bulk viscous coefficient in both 
of quark-gluon plasma and hadrons, which can be inserted in Eq. (|19p to 
calculate the expressions given in Eqs. (|20]l and (|2~T]l . 
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Fig. 3: At constant particle mass m and cross section a, rj is drawn againest 
T. There are two singularities at T — and at m/2. rj(m/T) is minimum at 
temperatures slightly below m/2. 



In the relativistic Boltzmann limit, the bulk viscosity \10\ [TT] reads 



which implies that rj doesn't depend on the extensive state quantity n, 
which gives the number of occupied states in momentum space, rj de- 
pends on the ratio of heat capacities and enthalpy per particle, which are 



T[K%(m/T) j 'K2(m/T)]/m, respectively. Also, it depends on mass m, tem- 
perature T and the cross section a. The latter has been given in Ref. |12j 
and is assumed to be constant for all states or particles. 

At m = 5 and constant <r, Eq. ([28]) is drawn in Fig. [3j It is clear that 
rj(m/T) has two singularities, one at T = and another one at m/2. It has 
a minimum value, at a temperature slightly below m/2. At much higher 
temperatures, rj(m/T) increases linearly with increasing T. In the high-T 
region, Eq. (|28p is likely no longer valid. 

In the Hagedorn picture, the particles, at very high energies, can be 
treated non-relativistically. In this limit, we express the partition function 
in mass spectrum p(m). 




(28) 



given by the auxiliary functions 





(29) 
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The mass spectrum p{m)dm give the number of states between m and m + 
dm. The non-relativistic part of the energy-momentum tensor is |13| 



(2vr)S 



pi pj pk pi 
mE 



-E/T 



p{m) dm d 3 p, 



(30) 



where r is the relaxation time. Taking into account the asymptotic behavior 
allows us to derive the bulk viscous term in T tJ , 



V = ^(2^ ln ( Cs_2 )' 



(31) 



where A is constant and c s = (dp/de) 1 ^ 2 is the speed of sound, which char- 
acterizes the propagating of signals in the cosmological background matter 
of the Early Universe. In the relativistic limit, the partial derivatives in c s 
are taken, adiabatically, i.e., at constant heat (or energy as we assumed in 
this model). 

In(c7 2 ) in Eq, @T) can roughly be estimated, when we approximate the 
thermodynamic quantities p and e, Eq. (l26j) . We assume that p and e are 
not changing with the bulk viscous coefficient 77, then 



m 




Using the dimensionless ratio m/T, last equation can be calculated, numer- 
ically. This is illustrated in Fig. [U The asymptotic value, ln(3), is fulfilled 
at high T. When T — > m, the function drops to a minimum value. It di- 
verges as long as T < m. The inverse of the relaxation time gives the drag 
coefficient of the background matter. In relativistic limit, we can model the 
relaxation time. For strongly coupled M = 4 SYM [TJ], r = /(1/T), 



ln(2) 



2vrT 



(32) 



To keep fitting scope of present work, we leave for a future work [9], the 
numerical estimates of Hubble paramter in the quantum cosmology with 
(section [5J and without (section [6]) bulk viscousity. 
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Fig. 4: At a constant particle mass, In (c s 2 ) is depicted as a function of T. 
Decreasing T lowers \n(c^ 2 ) below its asymptotic limit, ln(3). When T — > rn, 
In (c7 2 (m/T)) drops to a minimum value. 



7 Conclusions 

We have shown that the Universe, which is characterized by the FRW model, 
apparently obeys the laws of thermodynamics. We used classical assump- 
tions in order to derive the essential cosmological parameters, Hubble pa- 
rameter H, scaling factor a and curvature constant k. In doing so, we have 
assumed that the background matter is filled with an ideal thermal gas. Such 
a matter is homogeneously and isotropically distributed inside the available 
cosmological geometry. For simplicity, we assume that no interactions or 
phase transitions took place. 

The first gaol of this work is to study the effects of including finite bulk 
viscosity on the Early Universe using classical approaches. We started from 
the same assumptions as we did with the ideal thermal gas. We found con- 
siderable changes in all cosmological parameters. Comparing our results 
with the FRW model results in that the time-dependence of Hubble param- 
eter and scaling factor is slower than that of taking the background matter 
as an ideal thermal gas. 

Also, we have found a strong dependence of k on the thermodynamic 
quantities, like total energy E and bulk viscosity 77. The relation between E 
and the mass m determines the time scale, at which the negative curvature 
sets on. The total energy in the flat Universe is characterized by a dominant 
work of the bulk viscosity. When this takes place, the Universe decelerates. 
Otherwise, the expansion is positive. The expansion rate is directly propor- 
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tional to E and i]. As for thermodynamics of the Early Universe, we have 
found that the time evolution is affected by the bulk viscosity coefficient. 
Should this model be considered acceptable, essential modifications in the 
astrophysical observations are expected. 

The second goal of this work is to check the cosmological parameters in a 
quantum system. We started from basic assumptions of quantum mechanics 
and statistical physics. We expressed the Hubble parameter (and scaling 
factor in a straightforward way) in dependence on the infinitesimal changes 
in both phase and momentum spaces. We found that the Hubble parameter 
depends on 1/T and the time derivative of T. Based on this toy model, it is 
clear that the expansion of the Universe is derived by the generation of new 
states. In the relativistic limit, Hubble parameter depends on momentum 
space and the number of occupied state, besides the decay of T. 

Finally, we have studied the bulk viscosity in low-T and high-T regimes. 
For the first regime, the bulk viscosity diverges at vanishing T and at T ~ 
m/2, where m is mass. For the high-T regime, the bulk viscosity decreases 
with increasing T. Its asymptotic value is reached, when the speed of sound 
approaches its asymptotic limit. At vanishing T, the speed of sound diverges. 
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